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Abstract. To a unitary matrix U we associate a doubly stochastic matrix M by 
taking the modulus squared of each element of U. To study the connection between 
onset of quantum chaos on graphs and ergodicity of the underlying Markov chain, 
specified by M, we study the limiting distribution of the spectral gap of M when U 
is taken from the Circular Unitary Ensemble and the dimension N of U is taken to 
infinity. We prove that the limiting distribution is degenerate: the gap tends to its 
maximal value 1. The shape of the gap distribution for finite N is also discussed. 
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1. Introduction 

Recently it was proposed by Kottos and Smilansky 0] that quantum graphs can 
serve as a good "toy" model for quantum chaos. Numerical simulations confirmed that 
for some classes of graphs, including the fully connected graphs, the Bohigas-Giannoni- 
Schmit Conjecture || holds true. By this we mean that various statistical functions of 
the spectrum of the graphs converge to the corresponding functions obtained in Random 
Matrix Theory (RMT) [|J in the limit as the number of vertices is taken to infinity. 

On the graphs, the quantum evolution operator assumes the form of a finite unitary 
matrix U. In [|, @] it was argued that the classical counterpart of the quantum system is 
the Markov chain on the directed bonds of the graph with the matrix M of the transition 
probabilities given by 



It was conjectured that there is a link between ergodicity of the Markov chain and the 
convergence of the quantum statistics to RMT results. 

However, in |J it was shown that the statistics for a special family of graphs, the 
quantum star graphs, converge to a non-RMT limit as the size of the graphs increases. 
However "classically" the star graphs have nice properties: the corresponding Markov 
chain possesses a unique attracting equilibrium distribution which corresponds to the 
eigenvalue 1 of M. Thus ergodicity of each graph does not immediately imply RMT-like 
behaviour. 

It is well-known that a good quantitative measure of ergodicity of a Markov chain is 
its spectral gap. Given the spectrum {1, A2, A3, . . . Xn} of M, the spectral gap is defined 

as 



The spectral gap indicates the speed of convergence of any initial distribution to the 
equilibrium one: the greater g is, the faster is the convergence. Tanner || [7| suggested 
that if for a sequence of graphs the spectral gap of the corresponding Markov chains is 
uniformly bounded away from zero, then the spectral statistics will converge to their 
RMT form. This seems to be a good indicator of the conformance to RMT: while 
for each star graph the gap is non-zero, it decreases to zero as we increase the size of 
the graph. And in the numerically verified cases of RMT-like behaviour the gap stays 
bounded away from zero. 

To investigate the connection between random unitary matrices and the ergodicity 
of their "classical" analogues, we ask the following question: given a large random 
unitary matrix, find the probability distribution of the spectral gap of the corresponding 
Markov chain. 



Mb t v — \Ub t b' 



(1) 



= 1 — max IAJ. 

i=2,...N l 1 



(2) 
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2. Preliminaries 

Since the most natural ensemble of the unitary matrices is the Circular Unitary Ensemble 
(CUE), we shall restrict our attention to it, although the generalisation of the problem 
to other ensembles is straightforward. 

Definition 1. CUE(iV) is defined as the ensemble of all unitary N x N matrices 
endowed with the probability measure that is invariant under every automorphism 

U i — ► VUW, (3) 

where V and W are any two N x N unitary matrices. 

Definition 2. An entrywise nonnegative N x N matrix M is called doubly stochastic 
if 

N N 

]TM M = 1 Vj and E M m = 1 V *- ( 4 ) 

i=l j=l 

The set of all such N x N matrices we denote by DS(iV). 

If U is unitary then the matrix M defined by (|1|) is doubly stochastic. For a Markov 
chain it means that the uniform measure is invariant, that is (1/N, ... , 1/N) is a left 
eigenvector of M with the eigenvalue 1. 

Remark 1. It is easy to check that if A, B G DS(iV) then AB G DS(iV), thus DS(iV) 
forms a semigroup. 

Now we define the function S : CUE(iV) — > DS(iV) to map U to M according 
to ([[]). It is clear that S induces a probability measure on the semigroup DS(A r ). Thus 
our question is 

what is the probability distribution of the spectral gap g if the matrices M are 
selected randomly with the probability induced by the correspondence S. 

We denote the gap probability density function by f]y(g) and the corresponding 
cumulative distribution by Fj^(g) = Jq f^{g')dg'. For convenience we also denote by 
A2(M) the second largest (by modulus) eigenvalue of the doubly stochastic matrix M, 

|A 2 (M)| = .max |Ai| (5) 
and, as a consequence, g — 1 — |A 2 (M)|. 

3. Numerics and a conjecture 

First of all we are going to present the results of some numerical simulations. The 
algorithm we used is simple: random matrices from CUE(iV) are generated using 
the Hurwitz parametrisation || (see also 0), the spectrum of the corresponding 
doubly stochastic matrix M is computed and the gap is calculated according to the 
definition (|J). 
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The cumulative distribution functions obtained for N = 2, 3, 5, 10, 20, 40, 80 for 
sample sizes of 10 5 are shown on Figure [I]. From the plot we can see that the typical gap 
tends to 1 as we increase N. Thus it is natural to conjecture that the gap probability 
density function /jv(<?) converges to 8(g — 1) in the sense of distributions as N — > oo or, 
in other words, 

Pr {0 < g < a} -> as iV -> oo Vo < 1. (6) 

Further evidence in support of this guess is presented by Figure |2|, where the mean 
E| A2I and the standard deviation o of the second largest eigenvalue are plotted as 
functions of N. A good fit to the mean is given by E | A2 1 oc 1/yN (which is supported 
by heuristic observations presented in Section f|). Note that since < | A2 1 < 1 and 

E|A 2 | > aPr{a < |A 2 | < 1} + 0Pr{0 < |A 2 | < a} = a Pr {0 < g < a} , (7) 

to verify @ it is enough to show that E|A 2 | — > as N — > 00. 



4. Verification of the conjecture 

4.1. N = 2 case 

In the case of iV = 2 it is easy to calculate the gap probability explicitly. For 
completeness we include the derivation although it is obvious from Figure [I] that the 
distribution for N = 2 is not typical. A 2 x 2 unitary matrix can be parametrized as 

XJ = e iaf e^cos0 e ix sin0 \ 
I — e %x sin e~ 1 ^ cos <fi I ' 

with 

< a < 2tt, < x < 2tt, < tjj < 2vr, < <p < vr/2. (9) 

The Haar (uniform) measure on the CUE(2) can be expressed in terms of the parameters 
a, Xi an d 4> as (27r) -3 da d\ dip d sin 2 0. The corresponding doubly stochastic matrix 



is 



cos 2 d) sin 2 



M= ■ 2, 2\) ( 10 ) 

\ sm (j) cos z <f) J 

with eigenvalues 1 and cos 2 — sin 2 0. Thus the gap is equal to 2 min(l — sin 2 0, sin 2 0) 
and its distribution is uniform on the interval [0, 1] (which agrees with Figure |l|). 



4-2. Mean of the second largest eigenvalue 

Theorem 1. Let A 2 (M) be the second largest (by modulus) eigenvalue of the random 
doubly stochastic matrix M drawn from DS(iV) with the probability induced by the 
correspondence S : CUE(iV) — > DS(iV). Then 

E|A 2 |^0 as 00. (11) 

As a consequence, the gap distribution f^(g) converges to 5(g — 1) in the sense of 
distributions as N — > 00. 
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Proof. The eigenvalues of M G DS(iV) are, in general, complex. It would be more 
convenient to deal with nonnegative real eigenvalues. The trick is to consider the matrix 
A = M T M which, due to the semigroup property is also doubly stochastic. In addition, 
it is symmetric (all eigenvalues are real) and positive definite (i.e. all eigenvalues are 
nonnegative). To relate the eigenvalues of the matrices A and M we write for the second 
largest eigenvalue of A, 

A 2 (A) = max (Ax, a;) > (Ay,y) = (My, My) = |A 2 (M)| 2 , (12) 

x|=l, (x,e)=0 

where e = (1, . . . , 1) T is the eigenvector of both A and M with the eigenvalue 1, 
and y is the eigenvector of M corresponding to A 2 (M). From (|12|) and the inequality 
(E£) 2 < E(£ 2 ) we infer that to prove the Theorem it is enough to show that the mean 
of |A 2 (A)| converges to zero as iV — > oo. 

The simplest way to estimate the second largest eigenvalue of A is to compute its 
trace. Since A is positively defined, we have 

TrA n > 1 + A 2 (A) n . (13) 

Thus if we can find such n that E(Tr A n ) — > 1, it will imply that E(A 2 (A)) — > 0. It 
turns out that it is enough to take n = 2. However let us consider n = 1 too. 

N N N N 

Tr A = E(M T M) M = £ (M%(M),, = £ Mf. = £ \U hJ \\ (14) 
i=i *j=i *>i=i 2>i=i 

And calculating the mean 

N 

E(Tr A) = J2 El^l 4 = ^ 2 E|£/i,i| 4 , (15) 

where due to the invariance of measure of CUE(iV) the means of the different matrix 
elements are equal. To calculate E | tT" i , 1 1 4 one can either integrate over the measure of 



CUE(iV), as was done in ||10|| , or apply various invariance considerations II]. The result 
is E|[/i,i| 4 = 2/(N 2 + N) and therefore E(Tr A) -> 2 as N -> oo. This agrees with the 
numerical observation that E|A 2 (M)| oc 1/y/N: the eigenvalues of A are then of order 
1/N and there are N — 1 of them (not counting the 1). 
In the case n = 2 we have 

N N 

Tr A 2 = J2 (M T ) M -(M)^(M T ) M (M) ; , m = £ {U^U^U^ 2 

i,j,k,l=l i,j,k,l=l 

= ^2 \Uj,iUj,kUi,iUi t k\ 2 + ^2 Wj^u^if 

+ E l^,;M 4 + E l^l 8 - (16) 

ij^k,j=l i=k,j=l 

Applying the averaging we again find that due to the invariance of the measure all 
contributions from the first sum are the same and there are (iV 2 — iV) 2 of them; there 
are 2N 2 (N — 1) contributions from the second and the third sum (the contributions are 



Spectral gap of doubly stochastic matrices generated from CUE 



6 



equal) and N 2 contributions from the last. Counting the number of terms in each sum, 
we write 

ETr A 2 = (iV 2 - N) 2 E\U 1A U li2 U 2A U 2t2 \ 2 

+2N\N - l)E|Z7 1 ,iZ7 1)2 | 4 + iV 2 E|*7 M | 8 , (17) 



where the averages can be calculated using fllTf , 

E\u M vu»r N2 + N + 2 



E|t/i,it/i, 2 | 4 
E|f/i,i| 8 



N 2 (N 2 - l)(iV + 2)(JV + 3)' 
4 

N(N + 1)(N + 2)(N + 3)' 
24 



iV(iV + l)(iV + 2)(JV + 3)' 
Bringing everything together, we obtain 



ETr A 



N 3 + 8N 2 + UN - 2 „ / 1 



(N + 1)(N + 2)(N + 3) 
which effectively finishes the proof. 



(18) 
(19) 
(20) 

(21) 
□ 



5. Shape of the distribution of |A 2 (M)| 

It is now clear that the distributions /at ( 1 -^2 1 ) have singular limit. The natural question 
to ask, then, is what is the shape of /jv(|A 2 |) for finite, but large, N and whether there 
are sequences of normalizing coefficients and b^ such that the random variables 
(|A 2 | — a N )/b]y have continuous limiting distribution. While we are unable to give a 
definite answer, we can make a guess. Since |A 2 | = max 3= 2,...,Af |Aj|, it is not unnatural 
to conjecture that for large N the distribution of the second largest eigenvalue |A 2 | is well 
approximated by the Generalized Extreme Value distribution (GEV) [I2|, [T3|], specified 
by its cumulative distribution function 



G(x) 



exp < - 



X 



(22) 



for some (iV-dependent) £, a and b: the GEV distribution is used to model maxima of 
sequences of random variables. 

To check our conjecture we fit the numerical data for N = 80 to the form (^) using 
an Splus routine by Coles flT3"|. The routine minimizes the log-likelihood function 



l(a,b,£) 



In 6 



1+1/5 



X 



,(23) 



where x, is the sequence of observations of |A 2 |. The result of the routine is plotted as the 
probability density function G'(x) for fitted £ = —0.07, a and b against the histogram 
of the numerical data for |A 2 |. One can see that the agreement is very good. The 
probability and quantile plots, which we do not present here, also show good agreement. 
As N is taken to infinity, the values of a and b decay to zero, so that the limit of the 
distribution is the step function, and £ saturates at a nonzero value. 
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6. Conclusions 



We have shown that if we take element-wise modulus squared of a large unitary matrix 
U, the resulting doubly stochastic matrix M will have large spectral gap: with high 
probability all N — 1 "free" eigenvalues will have small modulus. In particular, it means 
that the Markov chain corresponding to M will quickly relax from any initial distribution 
into the uniform stationary state. 

While our observation is by no means a proof of Tanner's |7j conjecture that the 
families of graphs with large gap in their "classical" spectrum will adhere to RMT 
predictions, it is another powerful argument in its favour. We also note that the unitary 
evolution matrix of a general graph is sparse due to the topological restrictions [[l], ||]. 
Thus, although arbitrary U can be considered as the scattering matrix for an iV-star 
graph, the generalisation of the gap distribution to graphs with less trivial topology still 
remains to be studied. 

The results for the mean of the traces of M T M and (M T M) 2 , obtained in the 
course of our proof, completely agree with the numerical observation that the mean gap 
increases with N like 1 — (iV -1 / 2 ) . We also presented numerical evidence in favour 
of the claim that for finite N the distribution of the second largest eigenvalue is well 
approximated by the generalized extreme value distribution. 
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gap 

Figure 1. Cumulative distribution function of the gaps of the doubly stochastic 
matrices obtained from the unitary matrices from CUE(iV) for different values of N. 
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Figure 2. Estimations of the mean (circles) and the standard deviation (squares) of 
the modulus of the second largest eigenvalue A2(M) as functions of N. The error bars 
(based on the 95% confidence interval) are too small to be indicated: they are of order 
1.5% for the deviation and less than 0.2% for the mean. The lines correspond to the 
best fit with the power law. The exponents are —0.49 for the mean and —1.167 for the 
deviation. 
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Figure 3. Histogram of the numerical data from the simulation of | A2 1 for N = 80 
and the fitted GEV probability density function. 



